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Abstract: Polarized parton distributions and structure functions of the nucleon are ana- 
lyzed in the improved valon model. The valon representation provides a model to represent 
hadrons in terms of quarks, providing a unified description of bound state and scattering 
properties of hadrons. Polarized valon distributions are seen to play an important role 
in describing the spin dependence of parton distributions in the leading order (LO) and 
next-to-leading order (NLO) approximations. In the polarized convolution inte- 

gral is derived in the framework of the valon model. The Polarized valon distribution in 
a proton and the polarized parton distributions inside the valon are necessary to obtain 
the polarized parton distributions in a proton. Bernstein polynomial averages are used to 
extract the unknown parameters of the polarized valon distributions by fitting to the avail- 
able experimental data. The predictions for the NLO calculations of the polarized parton 
distributions and proton structure functions are compared with the LO approximation. It 
is shown that the results of the calculations for the proton structure function, xg^, and its 
first moment, T^, are in good agreement with the experimental data for a range of values 
of Q 2 . Finally the spin contribution of the valons to the proton is calculated. 
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1. Introduction 

Determination of parton distributions in a nucleoli in the framework of quantum chro- 
modynamics (QCD) always involves some model-dependent procedure. Instead of relying 
on mathematical simplicity as a guide, we take a viewpoint in which the physical picture 
of the nucleon structure is emphasized. That is, we consider the model for the nucleon 
which is compatible with the description of the bound state problem in terms of three 
constituent quarks. We adopt the view that these constituent quarks in the scattering 
problems should be regarded as the valence quark clusters rather than point-like objects. 
They have been referred to as valons. In the valon model, the proton consists of two "up" 
and one "down" valons. These valons thus carry the quantum numbers of the respec- 
tive valence quarks. Hwa jjj found evidence for the valons in the deep inelastic neutrino 
scattering data, suggested their existence and applied it to a variety of phenomena. Hwa 



- 1 - 



Q has also successfully formulated a treatment of the low-px reactions based on a struc- 
tural analysis of the valons. Some papers can be found in which the valon model has been 
used to extract new information for parton distributions and hadron structure functions Q . 

Hwa and Yang in two papers Q refined the idea of the valon model and extracted new 
results for the valon distributions. In a recent paper, unpolarized parton distributions and 
hadronic structure functions in the NLO approximation were extracted ||. The pion and 
meson structure functions in the valon model framework has also been analyzed || . In 
this paper we extend the idea of the valon model to both polarized and unpolarized cases 
in the LO and NLO approximations. 

The plan of the paper is as follows. In Subsec. 2.1, unpolarized valon distributions are 
explained. In Subsec. 2.2 we define a convolution integral, using a linear combination of 
spin- up and down quark distributions. In Subsec. 2.3 the relationship between unpolarized 
and polarized valon distributions in a proton is dealt with. The discussion of the polarized 
valons for the non-singlet and singlet cases is further refined in Subsec. 2.4. W functions 
which link the polarized valon distributions and the unpolarized ones are described in 
Subsec. 2.5. These contain some unknown parameters which can be obtained by fitting to 
the experimental data. In Sec. 3 the moments of the polarized valon and polarized parton 
distributions are analyzed in the NLO approximation. Subsec. 4.1 describes the method 
of Ref.|8| in which one constructs averages of the measured polarized structure functions 
weighted by suitably chosen Bernstein polynomials. These experimental Bernstein averages 
are then fitted in Subsec. 4.2, using the CERN subroutine MINUIT, to the QCD predictions 
for the corresponding linear combinations of moments. In section 5 the inverse Mellin 
transformation and the convolution integral are used to calculate the polarized parton 
distributions in the x-space. In Subsec. 6.1 the polarized structure function, xg\, for 
the proton is calculated in the LO and NLO approximations, and compared with the 
experimental data and some other theoretical models. Our prediction for the first moment 
of proton structure functions are presented in this subsection for some values of Q 2 . The 
spin contribution of valons to the proton are discussed and calculated in Subsec. 6.2. We 
also give our conclusions in Sec. 7. 

2. Valon distributions in unpolarized and polarized cases 
2.1 Valon distributions 

The subject of nucleon structure has mainly been investigated from two opposite view- 
points. On the one hand, one studies the bound state in terms of three quarks, such as 
in a bag model, and obtains various static properties of the nucleon and its spectroscopic 
partners. On the other hand, one probes the nucleon with high energy leptons in the hope 
of learning something about its three constituent quarks but finds structure functions that 
can only be understood in terms of an infinite number of partons (i.e., quarks, antiquarks 
and gluons). The two views can be reconciled if we recognize that the constituents in the 
quark model for the static nucleon are not the same objects as the quarks in the parton 



model. Since the scales of spatial resolution are different for the two problems, we may 
regard the former as clusters of the latter without contradicting either view. The clusters 
have the quantum numbers of the valence quarks, and have been called valons, for short. 
More precisely, a valon can be defined as a valence quark and associated sea quarks and glu- 
ons which arise in the dressing processes of QCD. In a bound state problem these processes 
are virtual and a good approximation for the problem is to consider a valon as an integral 
unit whose internal structure cannot be resolved. In a scattering situation, on the other 
hand, the virtual partons inside a valon can be excited and be put on the mass shell. It is 
therefore more appropriate to think of a valon as a cluster of partons with some momentum 
distribution. The proton, for example, has three valons which interact with each other in a 
way that is characterized by the valon wave function, while they respond independently in 
an inclusive hard collision with a Q 2 dependence that can be calculated in QCD at high Q 2 . 

In the static problems there is a little difference between the usual constituent quarks 
and the valons, since the point-like nature of the constituent quarks is not a crucial aspect 
of the description, and has been assumed mainly for simplicity. But, in the scattering 
problems it is important to recognize that the valons, being clusters of partons, can not 
easily undergo scattering as a whole. The fact that the bound-state problem of the nucleon 
can be well described by three constituent quarks implies that the spatial extensions of 
the valons do not overlap appreciably. A physical picture of the nucleon in terms of three 
valons is then quite analogous to the usual picture of the deuteron in terms of two nucleons. 

According to this description of the valon model, F^x^Q 2 ) as a structure function 
(e.g. F2 or X-F3) is related to a valon structure function corresponding to a valon v, by 
smearing of the valon momentum in the nucleon. On the other hand, the structure function 
of a hadron is obtained by convolution of two distributions: the valon distributions in the 
proton and the structure function for each valon. In the unpolarized situation we may 
write: 

FV(x,Q 2 ) = f dyG v/p (y)jq(-,Q 2 ) , (2.1) 

v J x y 

where the summation is over the three valons. Here G v / p {y) indicates the probability for 
the v- valon to have momentum fraction y in the proton. We shall also assume that the three 
valons carry all the momentum of the proton. This assumption is reasonable provided that 
the exchange of very soft gluons is responsible for the binding. Eq.(2.1) involves also the 
assumption that in the deep inelastic scattering at high Q 2 the valons are independently 
probed, since the shortness of interaction time makes it reasonable to ignore the response 
of the spectator valons. Thus, through Eq.(2.1) we have broken up the hadron structure 
problem into two parts. One part represented by G v / p (y), describes the wave functions 
of the proton in the valon representation. It contains all the hadronic complications due 
to the confinement. It is independent of Q 2 or the probe. The other part represented 
W ^(f^ 2 ): describes the virtual QCD processes of the gluon emissions and quark-pair 
creation. It refers to an individual valon independent of the other valons in the proton 
and consequently also independent of the confinement problem. It depends on Q 2 and the 



nature of the probe. 



Hwa and Zahir [0] assumed a simple form for the exclusive valon distribution which 
facilitated the phenomenological analysis as follows 

GuuD/ P (yuy2,y3) = giymTy^iyi + y 2 + y 3 - i) , (2.2) 

where yi is the momentum fraction of the i'th valon. The U and D type inclusive valon 
distributions can be obtained by double integration over the unspecified variables 



Gu/ P (y) = J dy 2 J dy 3 G UUD/p (y,y 2 ,y 3 ) (2.3) 
= gB(a + l,/3+l)y a (l-y) a+ P +1 , 

G D / P (y) = J ' dyi j dy 2 G UUD / p (y 1 ,y 2 ,y) (2.4) 
= gB{a + l,a + l) y P{l-y) 2a+l . 

The normalization parameter g has been fixed by 

/ G u/p (y)dy = [ G D/p (y)dy = 1 , (2.5) 
Jo Jo 

and is equal to g = [B(a + 1, f3+ l)B(a + 1, a + (3 + 2)] , where B(m, n) is the Euler-Beta 
function. 

R.C. Hwa and C. B. Yang [Q] have recalculated the unpolarized valon distribution in 
the proton with a new set of parameters. The new values of a , (5 are found to be a = 1.76 
and = 1.05. 



2.2 Convolution integral 

To describe the quark distribution q{x) in the valon model, one can try to relate the po- 
larized quark distribution functions or to the corresponding valon distributions 
and GK The polarized valon can still have the valence and sea quarks that are polarized 
in various directions, so long as the net polarization is that of the valon. When we have 
only one distribution q(x, Q 2 ) to analyze, it is sensible to use the convolution in the valon 
model to describe the proton structure in terms of the valons. In the case that we have 
two quantities, unpolarized and polarized distributions, there is a choice of which linear 
combination exhibits more physical content. Therefore, in our calculations we assume a 
linear combination of G^ and to determine respectively the unpolarized (G) and polar- 
ized (5G) valon distributions . 

To relate q^ and q^ to both G^ and G^, we can consider linear combinations as follows 



i /p ^Q 2 ) = E j\« G] /p (y) q lj/-,Q>) +P Gj /p(y )4(^Q 2 )]| , (2.6) 



and 



here and g^-l respectively denote the probability of finding q-up and (/-down in a G-up 
valon and etc. 

If we add and subtract the equations of (2.6, 2.7) we can determine the unpolarized 
and polarized quark distributions as follows 

Qi/p{x, Q 2 ) = qj /p (x, Q 2 ) + q\ /p {x, Q 2 ) 

= E jj G ] /P i< + + G ^h + Ptftfj-' ( 2 - 8 ) 

5q t/p (x,Q 2 ) = q} /p (x,Q 2 )-ql /p (x,Q 2 ) 

Since it is reasonable to assume that g^T = q^ and gT! = q^ , to obtain the unpolar- 
ized and polarized quark distributions in the proton, we need to choose in Eq.(2.8, 2.9) 
a = a' = = f3' = 1. 

Consequently we will have 

Q 2 ) = JjGj^Qi/j] + G^h/j])^, (2.10) 

and 

Sq i/P {x, Q 2 ) = E [\G] /p l^/ 3 ] ~ G l i/p [5q t/j })^-, (2.11) 

where q ifj = q}]. + q}^ = + q\j. and Sq ifj = qjj^ - qty. = qfo - are representing 
unpolarized and polarized quark distributions in a j- valon. Defining the unpolarized and 
polarized valon distributions as 

G 3tp - G] /p + G\ /p , (2.12) 



and 



we will obtain 



*G j/p = G] /p - G\ /v , (2.13) 



q i/p (x,Q 2 ) =^j\/j(^)GiMj > ( 2 - 14 ) 



and 



Sq i/p (x,Q 2 ) =J2 I S qi/j (-,Q 2 )6G j/p (y)^ . (2.15) 

■ J x y y 

As we can see similarly for the unpolarized case, the polarized quark distribution can be 
related to a polarized valon distribution. 

2.3 Spin dependence of the valon distributions 

The unpolarized and polarized quark distributions are respectively defined as 

q{x) = q (x) + q^(x) , 

5q(x) = q 1 * (x) — q^{x) . 

For the polarized parton distributions, \5f(x,Q 2 )\, and the unpolarized ones, f(x,Q 2 ), 
positivity requirements at low values of Q 2 imply the constraint [[l(], 11] 

\5f(x,Q 2 )\<f( x ,Q 2 ), (2.16) 

where / = u,u,d,d,s, s, g. Furthermore, we have the following sum rules [11] 



Au + An - Ad ~ Ad = A 3 = 1.2573 ± 0.0028 , (2.17) 

Au + Au + Ad + Ad - 2(As + As) = A 8 = 0.579 ± 0.025 , (2.18) 
AS = ^(Aq + Aq) = A 8 + 3(As + As) = A Q , (2.19) 



q 

where with n = 1, the first moment Af is defined by 



A/(Q 2 )= C dx5f(x,Q 2 ) . (2.20) 
J o 



It is helpful to consider the above reasonable theoretical constraints especially for 
the sea quark distributions. To determine the polarized parton distribution functions 
(PPDF's), the main step is to relate the polarized input densities to the unpolarized ones 
fll| 1 using some intuitive theoretical argument as the guideline. We employ the general 



ansatz for the polarized parton distributions of Ref.[ll], and introduce the following equa- 
tions to relate the polarized and unpolarized valon distributions 

SG j/p (y) = SW j (y)xG j/p (y), (2.21) 

here j refers to U and D type valons. The functions 5Wj(y) play an essential role in 
constructing the polarized valon distributions SGj/ p (y) from unpolarized ones, as we shall 
explain in Subsec. 2.5. 

Let us define the Merlin moments of any polarized U and D-valon distributions in the 
proton, SGu t D/ p (y), as follows: 

AM u/p (n) = f 1 y^SGu/piy) dy , (2.22) 
J o 

AM D/p (n) = f 1 y^SGn/^y) dy . (2.23) 
Jo 

Using polarized valon distributions, the moments of polarized valon distributions, 
AMu/ p (n) and AM D / p (n), can be easily calculated. The moments of polarized parton 
distributions in the proton can be written like the unpolarized case |J as follows 

Au v (n, Q 2 ) = 2AM NS (n, Q 2 ) x AM u/p (n) , (2.24) 

Ad v (n, Q 2 ) = AM NS (n, Q 2 ) x AM D/p (n) , (2.25) 
AE(n, Q 2 ) = AM s (n, Q 2 ) x (2AM u/p (n) + AM D/p {n)) . (2.26) 
The factor 2 in Eqs.(2.24, 2.26) is due to the existence of two-U type valons. 

In the above equations, AM NS (n, Q 2 ) and AM s (n, Q 2 ) are the moments of the non- 
singlet and singlet sectors which describe parton distributions inside the valon. We will 
introduce these functions in the LO and NLO approximations in Sec. 3. As we will see, in 
the LO approximation the first moment, AM NS (n = 1,Q 2 ) and AM s (n = 1,Q 2 ) will be 



equal to 1 for all values of Q 2 . So, using Eqs. (2.24-2.26) we have 

Au v = Au v (l, Q 2 ) = 2AM u/p {\) , (2.27) 

Ad v = Ad v (l, Q 2 ) = AM D/p (l) , (2.28) 

AS = AS(1, Q 2 ) = 2AM u/p {\) + AM D/p (l) . (2.29) 



It can be seen that in the LO approximation the first moment of the polarized quark dis- 
tributions is equal to the first moment of the related polarized valon distributions. We 
should notice that in the NLO approximation, the first moment of the polarized quark 
distributions are Q 2 -dependent, and the above equations are applicable only for Q 2 = Qq. 

The S symbol in Eq.(2.29) denotes Ylq=ud s(l + an< ^ consequently we have the 
following definition of AS 



A£(n,Q 2 ) = Y, (Aq(n,Q 2 )+Aq(n,Q 2 )), 

q=u,d,s 

therefore the moment of polarized sea quarks can be obtained as 



(2.30) 



Aq(n, Q 2 ) = (A£(n, Q 2 ) - Au v {n, Q 2 ) - Ad v {n, Q 2 ))/2f , (2.31) 

where / is the number of active quark flavours. If we choose n = 1 in Eq.(2.31) and 
use Eqs. (2.27-2.29), we can see that for all Q 2 values in the LO and for Q 2 = Qq in the 
NLO approximations, there is no contribution to the first moment of the sea quarks. As a 
result, we can not precisely determine the spin components of the hadrons. This difficulty 
is investigated in the following subsections. 



2.4 The improvement of polarized valons 

To resolve the above problem, first we improve the definition of polarized valon distribution 
functions 



.„ J SiWAy) x Gj/Jy) for non — singlet case , 
oGj/ p — > < . 



SWj (y) x Gj/ p (y) for singlet case 



(2.32) 



Using the above equation and according to the Eqs. (2.27-2.29), the first moment of the 
polarized u, d and £ distribution functions can be written as follows 

Au v = 2 / dy [SWu(y) x G u/p (y)} , (2.33) 
Jo 

Ad v = ^ dy [SW D (y) x G D/p (y)\ , (2.34) 
J o 



o 

DO 
O 
O 



O 

AS = J dy ([2 5W'u(y) x G u/p (y)] + [5W' D (y) x G u/p (y)]) . (2.35) 



The above equations can help us to consider the constraint of Eqs. (2. 17-2. 19) for the 
improved polarized valon model with an SU(3) flavour symmetry assumption. These con- 
straints play the same role as in the case of the unpolarized ones, Eq.(2.5), in controlling 
the parameter values which will appear in the polarized valon distributions. 

2.5 Description of the W function 

To define the actual y-dependence of the W' function, we parameterize this function as 

SW'jiy) = N jV a > (1 - yfi (1 + 7j y + r, jy - 5 ) . (2.36) 
As before the subscript j refers to U and i?-valons. 

The motivation for choosing this functional form is that the term y a i controls the low-y 
behavior valon densities, and (1 — y)P j that at large values y. The remaining polynomial 
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factor accounts for the additional medium- y values. 



For 5Wj (y) we choose the following form 

5 

SW-(y) = sWj{y) x Y, Amy^ ■ (2.37) 



m=0 



The additional term in the above equation, term), serves to control the behavior of 
the singlet sector at very low-y values in such a way that we can extract the sea quark 
contributions. Moreover, the functional form for 5W' and 5W" give us the best fitting \ 2 
value, as will be described in Sec. 4. In these functions, all of the parameters are unknown. 



Using the experimental data for g\ |L2|, 13] and the Bernstein polynomials, we can fit for 



the unknown parameters of Eqs.(2.36, 2.37). 

3. The NLO moments of PPDF's and the structure function 

In this section moments of the polarized parton distribution functions (PPDF's) in the 
valon are introduced. To calculate the Q 2 evolution of the polarized parton distributions, 
5f(x,Q 2 ), we need to Mellin transform these functions, i.e. 

Sf n (Q 2 )= C x n -Hf{x,Q 2 )dx . (3.1) 
J o 

The evolution of anomalous dimensions can be obtained by Mellin transformation of split- 
ting functions, and admit an expansion in powers of the running coupling constant a s {Q 2 ) 



Sd NS± = ^±lsA«h x ) + S^LI M W {l ) + ... , ( 3 .2 



whose detailed n-dependence has been specified in Ref.flO]. Here 5dqq is unique whereas 
all higher anomalous dimensions are scheme-dependent. 

The leading and next-to-leading order solutions of the renormalization group equation 
for the polarized moments can be expressed entirely in terms of the evolution parameter 



a s (Q 2 

The non-singlet (NS) moments evolve according to 



L{Q1 = • (3-3) 



AM^fn, Q 2 ) = (l - °" m ^°" iQo \ s4^ ~ L&dT ' (3-4) 



and the NLO running coupling constant is given by 



where b = 33 12 ^ an< ^ ^ = 15 |^ r 2 9 ■ in the above equations, we choose Qo = 1 GeV 2 as 
a fixed parameter and A is an unknown parameter which can be obtained by fitting to 
experimental data. 

The evolution in the flavor singlet and the gluon sector of the moments are governed 
by a 2 x 2 anomalous dimension matrix, with the following explicit solution 

AM s (n, Q 2 ) \ = f S i(°> , <Xs(Q 2 ) fT L 8dW« _ »s{Q 2 o) r5dW"fr\ (A ( o 6) 
AM 9q (n, Q 2 ) J V 2vr 2vr J I 1 J ' 1 ' ' 

where AM 99 indicates the spin dependent quark-to-gluon evolution function. All the as- 
sociated functions in the above equation have been defined in Ref.[14]. 

Using Eqs.(3.4, 3.6), the moments of NS, S and gq as a function of n are calculable. 
These moments in the LO and NLO approximations at Q 2 = 3 GeV 2 are presented in Fig. 
1. 

Determination of the moments of parton distributions in a proton can be done strictly 
through the moments of the polarized valon distributions. The parton distributions in the 
x-space which will be calculated, are 5u v , 5d v , 5T, and 5g. The moments of these distri- 
butions are denoted respectively by: Au v (n,Q 2 ), Ad v (n,Q 2 ), A£(n, Q 2 ) and Ag(n,Q 2 ). 
Therefore, the moments of the polarized u and d-valence quark in a proton are convolutions 
of two moments: 

Au v (n, Q 2 ) = 2AM NS (n, Q 2 ) x AM' u/p (n) , (3.7) 

Ad v (n, Q 2 ) = AM NS (n, Q 2 ) x AM D/p (n) , (3.8) 
The moment of the polarized singlet distribution (£) is as follows: 

A£(n,Q 2 ) = AM s (n,Q 2 )(2AMu /p (n) + AM' D/p {n)) . (3.9) 

For the leading and next-to-leading terms in Eqs.(3.4, 3.6), for the NS and S sectors, 
we can see that for all values of Q 2 in the LO and for Q 2 = Qq 2 in the NLO, the first 
moments of AM NS (1, Q 2 ) and AM S (1, Q 2 ) will be equal to one, and consequently the first 
moment of the polarized u and d quarks are respectively proportional to the first moment 
of the polarized U and D type valon distributions. 

For the gluon distribution we have 

Ag(n,Q 2 ) = AM^(n,Q 2 )(2AM' u/p (n) + AM D/p {n)) , (3.10) 
where AM 9q (n,Q 2 ) is the quark-to-gluon evolution function, given in Eq.(3.6). 



The NLO moment contribution of g%(x, Q 2 ) 



14] has the following form 



9i(n,Q 2 ) = + ^C^)[Aq(n,Q 2 ) + Aq(n,Q 2 )] + ^2^A 5 (n,Q 2 )} , (3.11) 

q 



where Aq(n,Q 2 ) = Aq v (n,Q 2 ) + Aq{n, Q 2 ), Aq(n,Q 2 ) and Ag(n,Q 2 ) are moments of 
the polarized parton distributions in a proton. <5C™, <5C™ are also the n-th moments of 
spin-dependent Wilson coefficients given by 



4 

bC n n = - 
q 3 



S 2 (n) + (S^n)) 2 - 1 5i(n) + 4r + 7T" + ' 



2 n(n + 1)7 ra 2 2n n + l 2 



(3.12) 



and 

^^[-^\)< s '< n » + 1 »-i + ^TT)l- < 3 ' 13 ' 

with Sk(n) defined as in Ref.[14]. 

So far the moments of the polarized parton distributions have been determined and 
the moments of the polarized proton structure function in the NLO approximation can 
be obtained by inserting the required distribution functions in Eq.(3.11). There are some 
unknown parameters in Eqs.(2.36, 2.37) and Eq.(3.5) which they also appear in Eqs.(3.7- 
3.10). It is obvious that the final form for g\(n,Q 2 ) involves the total of 17 unknown 
parameters. If the parameters can be obtained then the computation of all moments of the 
polarized parton distributions and the structure function, gi(n,Q 2 ), are possible. 



4. QCD analysis method 

4.1 The Bernstein averages of moments 

In the phenomenological investigations of the structure functions, for a given value of Q 2 , 
only a limited number of experimental points, covering a partial range of values of x, are 
available. Therefore, one cannot directly determine the moments. A method devised to 
deal with this situation is to take averages of the structure function weighted by suitable 
polynomials. We can compare theoretical predictions with the experimental results for the 
Bernstein averages, which are defined by |15| 



9n,k(Q 2 )= dxp nk (x)gi(x,Q 2 ) , (4.1) 
Jo 

where p n ,k{x) are the Bernstein polynomials, 

and are normalized to unity, f Q dxp n ,k{x) = 1. Therefore, the integral (4.1) represents an 
average of the function g\ (x, Q 2 ) in the region x n ^ — \ Ax n ±<x<x n ^ + \ Ax n ^ where x U: k 
is the average of x which is very near to the maximum of p n ,k( x ), an d Ax n fc is the spread 



of x n £. The key point is that values of g\ outside this interval have a small contribution 
to the above integral, as p n ,k{x) decreases to zero very quickly. In order to ensure the 
equivalence of the integral (4.1) to the same integral in the range x\ = x n ^ — ^Ax n ^ 
to X2 = x n> k + ^Ax n> k, we have to use the normalization factor of J^ 2 dxp n) k{x) in the 
denominator of equation (4.1) which obviously is not equal to 1. 

By a suitable choice of n, k we manage to adjust the region where the average is 
peaked to that in which we have experimental data |[^, 13 1. Using the binomial expansion 



in Eq.(4.2), it follows that the averages of g\ with p n ,k(x) as weight functions, can be 
obtained in terms of odd and even moments, 

and using the definition of Mellin moments of any hadron structure function we have 

*>* = rdVi^ + i) t SMk + + ^ ■ <44) 

We can only include a Bernstein average, g nj fc, if we have experimental points covering 
the whole range [x n ,k — \&.x n ^- l x n ^ + \Ax nt k] [16]. This means that with the available 



experimental data we can only use the following 41 averages, 



C P) (Q 2 ),C 2 xp) (Q 2 ), g£$m 

Another restriction which we assume here, is to ignore the effects of moments with high 
order n which do not strongly constrain the fits. To obtain these experimental averages 



from the E143 and SMC data [12, 13|, we fit xgi(x,Q ) for each bin in Q separately, to 



the convenient phenomenological expression 

xg [phen) = Ax B^ _ _ ^ 

This form ensures zero values for xg\ at x = 0, and x = 1. In Table I we have presented 
the numerical values of A, B and C for Q 2 = 3,5, 10 GeV 2 . A theoretical justification of 
Eq.(4.5) may be found in Ref.[17]. Using Eq.(4.5) with the fitted values of A, B,C one 



can then compute g^ p \Q 2 ) using Eq.(4.3), in terms of Gamma functions. Some sample 
experimental Bernstein averages are plotted in Fig. 2 in the LO and NLO approximations. 
The errors in the g^ p \Q 2 ) correspond to allowing the E143 and SMC data for xg\ to vary 
within the experimental error bars, including the experimental systematic and statistical 
errors [l3|] . We have only included data for Q 2 >3GeV 2 , this has the merit of simplifying 
the analysis by avoiding evolution through flavor thresholds. 



4.2 QCD fits to the Bernstein averages for gi 

Using Eq.(4.4), the 41 Bernstein averages g n ,k(Q 2 ) can be written in terms of odd and even 
moments g±(n,Q 2 ): 

g 2 ,i(Q 2 ) = 6{gi(2,Q 2 )- gi (3,Q 2 )) , 
52, 2 (Q 2 ) = 3( 5l (3,Q 2 )) , 

g^{Q 2 ) = 12 (51 (2, Q 2 ) - 24<7i (3, Q 2 ) + 12^(4, Q 2 )) , 
53 , 2 (Q 2 ) = 12( 5 i(3,Q 2 )-5i(4,Q 2 )) , 

5is,io(Q 2 ) = 4004.0 ( 5l (ll, Q 2 ) - ff i(14, Q 2 )) , 

-12012 ( 5l (12,Q 2 )- 5l (13,Q 2 )) . (4.6) 



The unknown parameters according to Eqs.(2.36, 2.37) will be, N(j,au,Pui • • • , ^o, ao 
Aq, Ai, ■ • • , Aq and finally A-^r from Eq.(3.5). Thus, there are 17 parameters to 
be simultaneously fitted to the experimental g n ,k(Q 2 ) averages. Using the CERN subrou- 
tine MINUIT []l8| , we defined a global \ 2 f° r an the experimental data points and found an 
acceptable fit with minimum x 2 /d.o.f. = 1.26 in the LO case and x 2 /d.o.f. = 0.94 in the 
NLO case with the standard error of order 10~ 3 . The best fit is indicated by some sample 
curves in Fig. 2. The minimum x 2 values for all 17 fitting parameters are listed in Table II. 

Now from Eqs.(2.36, 2.37), and using (2.32), we are able to determine the improved 
polarized U and D valon distributions in the proton. In Fig. 3 we have plotted the LO and 
NLO approximation results of ^ ■ y5Gj/ p for the non-singlet and singlet polarized valon 
distributions as a function of y. In this summation j runs over U, U, D valons. Since the 
polarized valon distributions are determined, we can obtain all polarized quark distributions 
in the proton using the convolution integral of Eq.(2.15). 

5. Polarized parton distributions in the x-space 

Polarized valon distributions in the proton have been calculated in Subsec. 2.4 and Sec. 
4. Using the polarized parton structure functions in a valon, it is then possible to extract 
polarized parton structure functions in a proton. To obtain the z-dependence of parton 
distributions from the n— dependent exact analytical solutions in the Mellin-moment space, 
one has to perform a numerical integral in order to invert the Mellin-transformation 

5f k ( z , Q 2 ) = - / dwlm[e^z- c - we AM k (n = c + we**, Q 2 )} , (5.1) 
^ Jo 

where the contour of the integration lies to the right of all singularities of AM^ (n = c + 
roe"^, Q 2 ) in the complex n-plane. For all practical purposes one may choose c ~ 1, (f> = 135° 
and an upper limit of integration, for any Q 2 , of about 5 + 10/ In instead of oo, which 



guarantees stable numerical results [19 



Inserting in Eq.(5.1) the three functions for k = NS, S, gq, i.e. AM NS (n,Q 2 ), 
AM s (n,Q 2 ) and AM 9q (n, Q 2 ) separately, we can obtain all polarized parton distribution 
functions inside the valon as a function of z, and at a fixed Q 2 value. We denote them by 
5f h (z = |,Q 2 ). These functions are independent of the type of valon but depend on Q 2 . 
As was discussed in Subsec. 2.2, the relationship between the polarized parton distributions 
of a proton and the polarized parton distributions in a valon can be given by Eq.(2.15) and 
now by using this equation, we can get the following expressions for the polarized parton 
distributions in a proton: 



5u v (x,Q 2 ) = 


2/V 5 (-,Q 
Jx y 


2 ) [sriuW 


x Gu/ P {y)\ — , 


5d v (x,Q 2 ) = 


/V 5 (-,Q 2 ) 
Jx y 


[5W' D (y) x 


G D / P (y)] j , 


5Z(x,Q 2 ) = 




2 [SW'u(y) 


xG u/p (y)] + [5W'n(y)xG D/p (y)]j 


Sg(x,Q 2 ) = 


/V*(-,Q 2 ) 
Jx y 


(2 [SW'uiy) 


xG u/p (y)] + [5W D (y)xG D/p (y)]) 



EE 

M 
hd 
l— 1 

dy 

y ' 
dy_ 
y 



(5.2) 



o 
o 

In Fig. 4 we have presented the polarized parton distributions in a proton at Q 2 = 3 GeV 2 . 
These distributions were calculated in the LO approximation and compared with some 
theoretical models [21-23]. In Fig. 5 we have presented the same distributions but in the 
NLO approximation. 



6. Results 



6.1 Polarized structure function 



In the last section PPDF's were calculated and now they can be used to extract the polar- 
ized structure function g\ in the LO and NLO approximations. 

According to the quark model, in the LO approximation, g\ can be written as a linear 



combination of 5q and 5q |lj, |24J , 



Q 2 ) = \Y. Q 2 ) + Sq(x, Q 2 )] , (6.1) 



2 



where e q are the electric charges of the (light) quark-flavours q = u,d, s. Here the sum usu- 
ally runs over the light quark-flavours q = u,d,s, since the heavy quark contributions (c, 
b,...) could preferably be calculated perturbatively from the intrinsic light quark (u,d,s) 
and gluon (g) partonic-constituents of the nucleon. 



o 

[\3 
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Within the MS factorization scheme the NLO contributions to g\{x,Q 2 ) are finally 



given by [14| 



+ ^p-[5C q ®{5q + 5q) + 25Cg®5g}} , (6.2) 



with the convolutions defined as 



(C®q)(x,Q 2 ) = [ ^-C(-)q(y,Q 2 ) . (6.3) 

J y y 

X 

The polarized proton structure function, xg\, in the LO and NLO approximations can 
now be presented. The results for xg\ in the LO and NLO approximations as a function of 
x and for some different values of Q 2 are depicted in Figs. 6,7. The comparison between dif- 
ferent models [21-23] and the available experimental data [12, 13 1 have also been presented. 

Since A(q + q) is the net number of right-handed quarks of flavor q inside a right- 
handed proton, it follows that ^AS is a measure of how much all quark flavours contribute 
to the spin of the proton. Similarly, Ag represents the total gluonic contribution to the 
spin of the nucleon. Using Eq.(2.20), first moments of the polarized parton distributions 
are calculated. These results have been shown in Table III. 

Finally we can define, similarly to Eq.(2.20), the first moment of g\ (the Ellis- Jaffe 
sum rule) by 

T\{Q 2 )^ f 1 dxg^Q 2 ) . (6.4) 
J o 

The results have also been given in Table III. 



6.2 Spin contribution of the valons to the proton 

The spin contribution of the valons to the proton can be analyzed through the valon model. 
The contribution of various polarized partons in a valon are calculable and by computing 
their first moment, the spin of the proton can be computed. In the framework of QCD the 
spin of the proton can be expressed in terms of the first moment of the total quark and 
gluon helicity distributions and their orbital angular momentum, i.e. 

\ = ±AY? + AgV + D> z , (6.5) 

where lT z is the total orbital angular momentum of all the quarks and gluons. 



Using Eq.(5.2) the contributions of 5T>(x,Q 2 ) and 5g(x,Q 2 ) in the proton can be 
calculated. These contributions for just one U valon in a proton can be extracted as 



8TF{x,Q 2 ) = f 6f s (*,Q 2 )[6W'u(y) x G u/p (y)]- , 

j x y y 

5g u (x,Q 2 ) = / 1 <5/"(-,Q 2 )[^(y) x G^y)]^ . (6.6) 

j x y y 

Using the above equations we can arrive at their first moments as in below 

AYy = f 5Z u (x,Q 2 )dx , 
Jo 

Ag u = [ 5g u (x,Q 2 )dx. (6.7) 
Jo 

The results for a D-valon can be obtained in a similar way. The resulting total quark 
and gluon helicity for a {7-valon is 

l -ATy + Ag u , (6.8) 

and for D-valon 

iAS D + Ag D . (6.9) 

Since each proton involves 2 [/-valons and one D-valon, the total quark and gluon helicity 
for the proton is 

^ATP + Af = 2(^AS f/ + A^) + ^AS D + Ag D . (6.10) 

With the assumption that the polarized U and D valons have spin of +^ and — ^ respec- 
tively, the determination of orbital angular momentum for each valon and finally for the 
proton, at different values of Q 2 , can be obtained from Eq.(6.5). Our results at Q 2 = Qq 
for the LO and NLO approximations are gathered in Table IV. 

7. Conclusions 

We have used the valon model to describe deep inelastic scattering. The model bridges the 
gap between the bound state problem and the scattering problem for hadrons. The valon 
distribution can serve as a model for the solution to the bound state problem. In contrast 
to this model, people usually use the GLAP equations to evolve the parton distributions 
from an initial value Qq. The valon model which was first introduced by Hwa [1-7,9], gives 
us a clear insight as to how to construct the hadron structure from the parton distributions. 
According to this model, a valon is defined as a cluster of valence quarks accompanied by a 
cloud of sea quarks and gluons. It can be considered as a bound state in which for instance 
a proton consists of three valons, two U-valons and one D-valon. This model has been 
applied to the unpolarized case and gives acceptable phenomenological results. 



In this paper we extended the idea of the valon model to the polarized case to describe 
the spin dependence of the hadron structure functions. The polarized valon distribution is 
derived from the unpolarized valon distribution. Prom the phenomenological point view, 
we need to define weight functions to construct a functional form for the polarized valon 
distribution. Due to the difficulty which exists in deriving the spin contribution of the 
sea quarks for all values Q 2 in the LO and at Q 2 = Qq in the NLO approximations, we 
improved the weight functions and introduced two different types of them for the singlet 
and non-singlet cases. In deriving the polarized valon distributions some unknown param- 
eters are present which should be determined by fitting to the experimental data. Here we 
have used a method similar to Ref.[||, to fit QCD predictions for the moments of the g± 
structure function to the suitably constructed Bernstein polynomial averages of the E143 
and SMC experimental data. 

After calculating polarized valon distributions and all the parton distributions in a 
valon, we calculated the polarized parton density in a proton. The results were used to 
evaluate the spin components of the proton. It turns out that our results for the polarized 
structure functions are in good agreement with all the available experimental data of g\ 
for the proton. Our prediction for I^, the first moment of g\, is also in good agreement 
with the experimental data. 

We hope to report in future on application of the improved valon model in describing 
more complicated hadron structure functions. We also hope to be able to consider the 
symmetry breaking of polarized sea quarks in our new calculations. 
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Table Caption 



Q 2 GeV 2 


A 


B 


C 


3 


0.274 


0.876 


1.212 


5 


0.339 


0.911 


1.724 


10 


0.342 


0.842 


1.932 



Table I Numerical values of fitting A, B,C parameters in Eq.(4.5). 



LO 



j 


N 3 


(Xj 




7i 


Vj 




u 


0.002 


-2.3789 


-1.7518 


11.0804 


-1.4629 




D 


-0.005 


-1.5465 


-1.8776 


8.5042 


-0.8608 






A 


A l 


A 2 


^3 


A 4 


A 5 




0.00043 


0.2954 


-6.9134 


30.9851 


-39.7383 


16.4605 


203 



NLO 



j 




dj 


Pi 


7i 


Vj 




u 


0.0038 


-2.1501 


-0.8859 


10.6537 


-0.1548 




D 


-0.0046 


-1.5859 


-1.5835 


9.6205 


-0.8410 






A 


Ai 


A 2 


^3 


M 


A 5 


A MS ( MeV ) 


-0.0025 


-3.1148 


15.8114 


-21.1500 


10.5025 


-0.9162 


235 



Table II Numerical values of fitting parameters for the best fit of Fig. 2 in the 
and NLO approximations. 



Q 2 (GeV 2 ) Au v Ad v AS Aq Ag 

1 0.8807 -0.3328 0.1786 -0.0615 0.5480 0.1214 

3 0.8774 -0.3315 0.1724 -0.0622 0.7774 0.1242 

5 0.8765 -0.3311 0.1706 -0.0624 0.8737 0.1249 

10 0.8754 -0.33308 0.1688 -0.0626 0.9986 0.1258 



Table III The first moments of polarized parton distributions, Au v , Ad v , AS, Aq, 
and r i in the NLO approximation for some value of Q 2 . 
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IAS* 


Ag* 


K 


iAS i + Ag i + L\ 


U 


0.091 


0.440 


-0.031 


0.5 


D 


-0.077 


-0.332 


-0.091 


-0.5 


P 
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0.548 


-0.153 


0.5 



(a) 



i 


IA^ 


Ag* 


4 


±A£* + Ag l + L\ 


U 


0.039 


0.440 


0.020 


0.5 


D 


0.011 


-0.332 


-0.178 


-0.5 


P 


0.089 


0.548 


-0.137 


0.5 



(b) 

Table IV Spin contribution of valons and proton in the LO (a) and the NLO (b) 



approximations at Q = Q 



o- 
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Figure Caption 



Fig. 1-The contribution of moments NS, S and gq from Eqs.(3.4, 3.6) as a function 
of n at Q 2 = 3 GeV 2 . Dashed-dotted line and solid line indicate respectively the LO and 
the NLO approximations. 

Fig. 2-Fit to xg\ using the Bernstein averages. Dashed-dotted line and solid line in- 
dicate respectively the LO and the NLO approximations. 

Fig. 3- The plots of total polarized valon distributions in the proton, ^ • y5Gj/ p , for 
non-singlet and singlet polarized valon distributions as a function of y where j runs over 
U, U, D valons. 

Fig. 4-Polarized parton distributions in the proton at Q 2 =3 GeV 2 as a function of x 
in the LO approximation. The solid line is our model, dashed line is the AAC model 
(ISET=1)[21], dashed-dotted line is the BB model (ISET=1)[22] and long-dashed line is 
the GRSV model (ISET=3)[23]. 

Fig. 5-Polarized parton distributions in the proton at Q 2 =3 GeV 2 as a function of x 
in the NLO approximation. The solid line is our model, dashed line is the AAC model 
(ISET=3)[21], dashed-dotted line is the BB model (ISET=3)[22] and long-dashed line is 
the GRSV model (ISET=1)[23]. 

Fig. 6-Polarized proton structure function xg\ as a function of x which is compared with 
the experimental data from Ref. [12,13] for different Q 2 values in the LO approximation. 
The solid line is our model, dashed line is the AAC model (ISET=1)[21], dashed-dotted 
line is the BB model (ISET=1)[22] and dashed-dotted-dotted line is the GRSV model 
(ISET=3)[23]. 

Fig. 7-Polarized proton structure function xg\ as a function of x which is compared with 
the experimental data from Ref. [12,13] for different Q 2 values in the NLO approximation. 
The solid line is our model, dashed line is the AAC model (ISET=3)[21], dashed-dotted 
line is the BB model (ISET=3)[22] and dashed-dotted-dotted line is the GRSV model 
(ISET=1)[23]. 
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